Generalized virial theorem in Palatini f{1Z) gravity 



A. S. Sefiedgar* K. Atazadeh''' and H. R. Sepangi ^ 
Department of Physics, Shahid Beheshti university, Evin, Tehran 19839, Iran 

August 26, 2009 

'g< 

O 
O 

(N 

.' Abstract 

^ , We use the colhsion-free Boltzmann equation in Palatini f{TV) gravity to derive the virial the- 

' orem within the context of the Palatini approach. It is shown that the virial mass is proportional 

, to certain geometrical terms appearing in the Einstein field equations which contribute to grav- 

^SJ ■ itational energy and that such geometric mass can be attributed to the virial mass discrepancy 

in cluster of galaxies. We then derive the velocity dispersion relation for clusters followed by the 
^ ■ metric tensor components inside the cluster as well as the f{TV) lagrangian in terms of the obser- 

vational parameters. Since these quantities may also be obtained experimentally, the f{TZ) virial 
^ • theorem is a convenient tool to test the viability of f{TZ) theories in different models. Finally, 

{3J[)| we discuss the limitations of our approach in the light of the cosmological averaging used and 

questions that have been raised in the literature against such averaging procedures in the context 
I of the present work. 
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Dark matter is presently one of the most exciting open problems in cosmology. There are some com- 
pelling observational evidence for the existence of dark matter for which, the galaxy rotation curves 
and mass discrepancy in cluster of galaxies are two prominent examples. According to Newtonian 



\ gravity, galaxy rotation curves give the velocity of matter rotating in a spiral disk as a function of 

the distance from the center of galaxy according to v{r) = GM{r) jr. If we assume that the cluster 
mass obeys the relation M(r) = p^^^, where p is considered as a constant density in the cluster, 
then the velocity increases linearly within the cluster and drops off as the square root of r outside 
the cluster. However, observation shows that the velocity remains approximately constant, that is 
M ~ r. This points to the possible existence of a new invisible matter which is referred to as dark 
matter and is distributed spherically around galaxies [1]. 

What is now known as the mass discrepancy of clusters can be understood when estimating the 
total mass of a cluster in two different way; cluster masses can be deduced by summing the individual 
member masses which we shall call M in total. Alternatively, the virial theorem can be used to 
estimate the mass of a cluster, My, by studying the motions of each member of the cluster. As it 
turns out. My is nearly 20 — 30 times greater than M and this difference is known as the virial mass 
discrepancy [1]. The prime tool in dealing with the above discrepancy is to postulate dark matter. 
There are several candidates for dark matter. One possible categorization tells us that dark matter 



*e-mail: a-sefidgar@sbu.ac.ir 
^e-mail: k-atazadeh@sbu.ac.ir 
*e-mail: hr-sepangi@sbu.ac.ir 



1 



can be baryonic or non-baryonic. The main baryonic candidates are the Massive Compact Halo 
Objects (MACHO) which include brown dwarf stars and black holes. The non-baryonic candidates 
are basically elementary particles which have non-standard properties. Among the non-baryonic 
candidates, we can point to axions as a solution to the strong CP problem. However, the largest 
class is the Weakly Interacting Massive Particle (WIMP) class which consists of hundreds of, as yet, 
unknown particles [2]. The most popular of these WIMPs is the neutralino from supersymmetry. 
WIMP's interaction cross-section with normal baryonic matter is extremely small but non-zero, so 
their direct detection is a possibility. Neutrinos, may also be considered as possible candidates for 
dark matter. Another important categorization tells us that dark matter can be hot or cold. A dark 
matter candidate is called hot if it was moving at relativistic speeds at the time when galaxies were 
just starting to form. It is called cold if it was moving non-relativistically at that time. Of the above 
candidates only the light neutrinos would be hot, all the others would be cold. There is, as of now, 
no non-gravitational evidence for dark matter. Moreover, accelerator and reactor experiments do not 
support the scenarios in which dark matter emerges. 

To deal with the question of dark matter, a great number of efforts has been concentrated on 
various modifications to the Einstein field equations [3]. One such modification is that of f{TZ) where 
TZ is the Ricci scalar. Theories of f(JZ) modified gravity have had some success in explaining the 
accelerated expansion of the universe [4, 5] and account for the existence of dark matter [6, 7, 8]. In 
this paper we study f{Tl) gravity in the context of the Palatini formalism. As is well known, starting 
with the usual Einstein-Hilbert action, both the Palatini and metric approaches result in the same field 
equations. However, if the action is taken as a generic function of TZ, then the two approaches result in 
different field equations [9]. Here, we study the virial theorem within the framework described above. 
In general, virial theorem plays an important role in astrophysical objects like galaxies, clusters and 
super clusters. By using the virial theorem and studying the observational data from the velocity 
of each member, one can estimate the mean density of such objects, rendering the prediction of the 
total mass possible. The virial theorem also offers interesting predictions on the stability of the 
astrophysical objects. Several authors have studied the virial theorem in models with a cosmological 
constant [10, 11], branc world scenarios [12] and metric /{TZ) theories [13]. Our main purpose in this 
paper is to obtain the generalized form of the virial theorem in f(JZ) Palatini formalism by using the 
collisionless Boltzmann equation. Of course, some extra terms emerge in virial theorem which are 
originated from the modified action and are geometric in nature. We show that one may account for 
the virial mass discrepancy by taking into account such extra terms. The components of the metric 
tensor inside the galaxies may also be derived in terms of physical observable quantities like the 
temperature of intra cluster gas and radius and density of the cluster core. Finding the components 
of the metric tensor leads to the form of the Lagrangian in f{TZ) gravity in terms of observable 
quantities. Thus, the virial theorem is a convenient tool to test the validity of f(JZ) models. 

In what follows, we first give a brief review of the Palatini formalism and the gravitational field 
equations are derived in scalar tensor representation of f{TZ) gravity. Next, we introduce the relativis- 
tic Boltzmann equation from which we deduce the virial theorem with the aid of the field equations 
and discuss the limitations of the assumptions made. The geometric mass and density of a cluster is 
then identified in terms of the observable quantities and the metric components are calculated inside 
the cluster. We then move on to study the velocity dispersion relation in galaxies. Finally, we present 
the Lagrangian which represents the f{TZ) theory inside the cluster. Conclusions are drawn in the 
last section. 

2 Palatini f{7l) gravity in scalar-tensor representation 

Let us start with the action 

Spalatini = [ d'^X^/^f{^Z) + f d!^Xy/^jCrnigtMv), (1) 
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where is the matter lagrangian and k = SttG. The Ricci scalar is written as TZ in the context 
of the Palatini formahsm to point out that it is different from the Ricci scalar, R, in the context of 
metric f{R) gravity. It is necessary to stress that Tl = g^^'^TlnviV), where 7?.^jy(r) is constructed from 
the connection which is independent of the metric. In addition, in Palatini approach, the lagrangian 
corresponding to matter does not depend on the connection. Varying the action with respect to the 
metric and connection, respectively, yields 

F(7^)7^^, - ^f{n)g^, = -kT^,, (2) 



-V^{V^F{n)gn + VAV^Frng'^^nS"^ = 0, (3) 

where we have used STZ^i, = Va^P^^ — Vj^^F^;^ with being the covariant derivative which is defined 
with the independent connection and (fiv) and [/xi/] define the symmetric and anti-symmetric parts 

of the relevant parameter and T^^ is the energy-momentum tensor. We also denote F(TZ) = ^/(^)- 
By contracting equations (2) and (3), the field equations can be written as 

F(7^)7^ - 2/(7^) = -kT, (4) 



-Vx{V^F{n)g^n = 0- (5) 
It is useful to define a metric conformal to gfj,,^ as follows 

V = Fin)g^„ (6) 

which yields 

V^h'^'' = ^gFin)g''r (7) 

Now, equation (3) becomes the definition of the Levi Civita connection of hfj^n and can be solved 
algebraically to give 

r^. = ^h}" [d^^Ka + 5. V - daKA ■ (8) 

By using equation (6), we can write the connection in terms of g^^ 

= + MF{n)g^a) - d,{F{n)g^,)\. (9) 

We now have an expression for P^j, in terms of TZ and ^^jy so the independent connection can be 

eliminated from the field equations. In fact, using conformal transformation in equation (6), the 
Riemann tensor, Ricci scalar and Einstein tensor in Palatini formalism can be derived in terms of the 
metric ones [14] 

7^M. = R^.. - ^ + -^)^^^^^'^ + \9^..^)^{n). (10) 

Contracting equation (10) by g^^ yields the Ricci scalar 

n = R- \^,v,Fin)vmn) + ^□F(7^). (ii) 

Substituting equations (10) and (11) in (2), we may calculate the Einstein tensor 



G^u — TpFfiv + 
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As a result, the independent connections disappear and the theory is brought to the form of GR with 
a modified source which only depends on the metric and matter fields. We are now ready to introduce 
a Legendre transformation {TZ, /} {(f), V} defined as 

= F(7^), y(0)^7^(0)F-/(7^(</.)). (13) 

The theory in the new representation is given by the action 

Spalatini = J d!^Xy/^[(t)'R - 1/(0)] + j (f Xy/^Cm{9 ij.u) ■ (14) 

Let us write the field equations in terms of the new parameters [9] 
where 

:iy^.(t>^u(^ - ^5/.^VVVa<^) - ^(V;.V^ - g^yU)cl) - ^g^^. (16) 



It is now easy to realize that this result is the same as that of the Brans-Dicke theory with a; = |. 

Introducing a modified action, f{TZ), results in the appearance of an effective gravitational constant. 
Note that Keff = ^ = ^'^Geff- Therefore, we may find from equation (15) that Geff = ^. The 6^^ 
tensor on the right hand side of equation (15) emerges as a new additional source for the gravitational 
field. 



3 Field equations for a system of identical and collision-less point 
particles 

Let us now consider an isolated and spherically symmetric cluster being described by a static and 
spherically symmetric metric 

= -e''^')dt^ + e^(")dr2 + r'^dO'^ + sin^ Qd^p^ . (17) 

Suppose that the clusters are constructed from identical and collision-less point particles (galaxies), 
being described by the distribution function /b. The energy-momentum tensor may be written in 

terms of /s as [15] 

Tjxv = J fsmUfj^Ui^du, (18) 

where m is the cluster's member mass, u is the four velocity of the galaxy and du = '^"'•'^"^^"y jg 
the invariant volume element of the velocity space. The energy momentum tensor of the matter in a 
cluster can be represented in terms of an effective density Pejf and an effective anisotropic pressure, 
with radial p^^j and tangential p-^j^j^ components [10]. In other words, we have 

Peff = p{ut), ptff = pWr), piff = p{ue) = pWI), (19) 

where ( ) represents the usual macroscopic averaging. 

Before going any further, a word of caution is in order at this point. There is an interesting 
subtlety concerning the energy momentum tensor appearing on the right hand side of the Einstein 
field equations defined in terms of the averaged velocities; the fact that one can average over ve- 
locities to obtain an averaged energy-momentum tensor docs not necessarily mean that there exists 
a corresponding averaged space-time metric associated to that averaged energy-momentum tensor. 
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This point was first raised by Flanagan [16] and Olmo [17]. Later on, Motta and Shaw, using genal 
averaging arguments [18] claimed that the averaged metric is almost indistinguishable from that of 
General Relativity. Subsequently however, it was shown in [19] that such conclusions are not always 
correct because of the existence of counter examples; infrared corrected models do not allow such 
averaging of the metric. Models with high energy corrections do admit such averaging, but it is not 
guaranteed in general. 

The above arguments boil down to the fact that in modified gravity theories such as f{Tl) and 

within the context of the Palatini approach, it may not be reasonable to replace the microscopic 
energy-momentum tensor by its cosmological average. Modified gravity theories imply non-linear 
correction terms to the energy-momentum tensor. When applying the field equations to cosmological 
scales, one may consider the microscopic structure of all particles in the system. The existence of 
non-linear terms means that one should average over all the microscopic structure of matter and 
it seems as if in cosmological scales the usual macroscopic averaging procedure is no longer valid. 
In Einstein general relativity, the field equations are approximately linear on microscopic scales. 
Therefore, the microscopic structure of matter is not particularly important on macroscopic scales. 
The condition in metric f{R) gravity is the same as in Einstein general relativity. But what about the 
Palatini formalism? The Palatini formalism has been applied to f{TZ) gravity in the Einstein frame 
[16, 17, 19] and the field equations thus obtained are non-linear in terms of the energy-momentum 
tensor even on the smallest scales. In fact, this non-linearity leads to the problem of averaging in 
Palatini theories. In other words, the standard averaging procedure may no longer be valid and results 
in incorrect predictions. In addition, if one takes into account all the microscopic structure of the 
matter, the palatini formalism of f(TZ) gravity will be indistinguishable from the standard General 
Relativity with a cosmological constant [18]. Interestingly, although Palatini theories were designed to 
modify gravity on large scales, they actually modify physics on the smallest scales and leave the large 
scales practically unaltered. In this work, we study clusters containing collision-less point particles 
(galaxies) within the Palatini formalism in the Jordan frame. Here we take the cosmological averaging 
of the energy momentum tensor without studying the microscopic structure of the matter. Of course, 
the validity of our results should be taken in the light of the discussion presented above. 

Using r^jy = {p + p)unUi, + pg^i, and u^U/^ = — 1 with above definitions, the field equations become 



+^(V*V.-n)<^+^, (20) 



7 + ^:2 



4(vv, 



2VVVa</)) 



(21) 



'- V u'X' v" 



_l(V^V.-D)^-i^, 



(22) 



v' - A' v'y v" 



2r 



4 +Y + T 



1 , V 



^(V-V,-n)^--. 



(23) 
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Another useful equation is obtained by summing the above four equations 



2v' v'>! „ v'^ 

Vv ^ 

r 2 2 



~p{u^) - ^(2Vi0VV) + ^(2V*V* + 0)0 - ^, (24) 



where («^) = (u^) + (w^) + (n^) + (m^). Since we are interested in the extra-galactic region, we assume 
a small deviation from standard general relativity in which we have ^=\. Let us take ^ = 1 + eg'{Tl) 
in our case where e is a small quantity and g'(JZ) describes the modification of the geometry due to 
the presence of tensor 0^^^ [20] . Using 1/^ ~ ^ — eg' ij^) , equation (24) can be written as follows 



2iy' v'\' ,, z/ 



,/2 



-Kp{v?) - Kp^, 



where 



It is not difficult to see that equation (26) can be written as 

,/2" 



=l+eg'{TV) 



(25) 



(26) 



u' ly'X' u" V 



-47rGp(«^) - 47rGp<^. 



(27) 



4 The virial theorem in Palatini /(7^) gravity 

To derive the virial theorem, we need the Boltzmann equation which governs the evolution of the 
distribution function. Integrating this equation on the velocity space when accompanied by the grav- 
itational field equations can yield the virial theorem. We consider an isolated spherically symmetric 
cluster which is described by equation (17). The galaxies in the cluster behave like identical, collision- 
less point particles. The distribution function is denoted by Jb which obeys the general relativistic 
Boltzmann equation. Space-time is a time oriented Lorantzian four dimensional manifold. The tan- 
gent bundle T{M) is a real vector bundle whose fibers at a point x G M is given by the tangent space 
Tx{M). The state of a particle is given by the four momentum p G T!r(M) at an event x G M. The 
one particle phase space Pphase is a subset of the tangent bundle given by [13, 15] 



phase •" 



{{x,p)\x G M,p G T^{M),p^ = -ml), 



(28) 



where tuq is the particle mass. A state of a multi particle system can be described by a continuous 
non-negative function f{x,p). It is defined on Pphase and gives the number dN of the particles 
crossing the volume dV with momenta p lying within a corresponding three-surface element dp in the 
momentum space. The mean value of /b equals to the average number of occTipicd particle states 
{x,p) [13, 15]. Let {x"} with a = • • • 3 be a local coordinate system in Af , defined in some open 
set U C M. Note that dt is timelike future directed and da, a = 1,2,3 are spacelikc. Then {gf^r} is 



p"m- We 
A 



the natural basis for tangent vectors. Each tangent vector in U can be written as p 
can define a system of local coordinates {z^}, A = 0, 7 in Tu{M) as = a;" and 2"+^ = p" 
vertical vector field over T{M) is given by tt = p"'-^- The geodesic field cr, which can be constructed 
over the tangent bundle, is defined as 

^-P P P i^a^Q^p -P i^a, 

where Pa-y are the connection coefficients. Physically, a describes the phase fiow for a stream of 
particles whose motion through spacetime is geodesic. Therefore the transport equation for the 
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propagation of a particle in a curved arbitrary Riemannian spacetime is given by the Boltzmann 
equation [13, 15] 

where i = 1,2,3. In many apphcations, it is convenient to introduce an appropriate orthonormal 
frame or tetrad e^{x), a = - ■ - 3 which varies smoothly over some coordinates in the neighborhood 
of U and satisfies the condition g^^e'^^e\, = ij"'^ for all x e U. Any tangent vector p'* at x can be 
expressed as p'^ = p"ea, which defines the tetrad components p". In the case of the spherically 
symmetric line element given by equation (17), we introduce the following frame of orthonormal 
vectors [13, 15] 

e°=e'^/2^o, el = e^%, el = rS^ el = Tsiue5% (30) 

where u*^ is the four velocity of a typical galaxy, satisfying the condition u^^u^ = — 1 with tetrad 
components = u'^e^. The relativistic boltzmann equation in tetrad components can be written as 

''<0^ + ^^ouu^ = 0, (31) 

where the distribution function /b = fB{x'^,u'^) and 7^^ = e'^.^e'^e'^ are the Ricci rotation coefficients 
[10, 13, 15]. We may assume that /b depends only on the radial coordinate r. Using equation (9), 
the relativistic Boltzmann equation in Palatini formalism is obtained as 

„ . dfB _ (l,,2dv _ '4+^1 
"1 dui \2^0dr r 



EL 
2F 



It - + -3|£) - ^-3 cot ^eV2 _ ,3|g) 

[ul + ul-ul-ul)^ + 2u,U2^ + 2nin3|gj = 0. (32) 



Since we have assumed the system to be spherically symmetric, the term proportional to cot 9 must 
be zero. Let us take [15] 

UO = Ut, Ui = Ur, U2 = U0, U3 = U^. (33) 

Multiplying equation (32) by murdu, integrating over the velocity space and assuming that /b van- 
ishes sufficiently rapidly as the velocities tend to ±00, we obtain 

+§[(u?)+9(u^)-(u^)-«)]=0. (34) 

It is worth mentioning that in integrating equation (32) we note that the functions F and F' depend 
on the average of the square of velocities according to equation (19) via equation (4), which are 
assumed to be constant quantities. Now, it is useful to multiply equation (34) by 47rr^ and integrate 
over the cluster volume to obtain 



I 



£ [(^?) + 9(^?) - i'^e) - K)] ^^r'dr = 0, (35) 



where R is the radius of the cluster. 
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At this point, it is appropriate to introduce some approximations. First, consider that A' and v' 
are small quantities. Then the terms proportional to u'X' and v'"^ in equation (27) may be ignored. 
Then, assuming that e"^ 1 inside the cluster [13], we can write equation (27) as 

2^l^("'l^)="^^^^"^^^^'^- ^^^^ 



Second, consider that the galaxies in the cluster have velocities much smaller than the velocity of 

4) 



light. In other words, («^) {uq) (n^) ^ Ri 1. Now equation (35) can be written as 



1 du F' 

2K + - p—Anr^dr + / —pATrr^dr = 0, (37) 

2 Jo c'r Jo 



where 

rR 

K = - + {uj) + «)]27rr2dr, (38) 

Jo 

is the kinetic energy of the galaxies. Multiplying equation (36) by and integrating yields 

GM{r) = "2^'^ - GM^ir), (39) 

where we have used M = dM{r) = Airpr^dr as the total mass and we have also defined 
M(f,{r) = Att Jq P4>{r')r''^dr' as the geometric mass of the system. Now, consider the definitions 

n = -r^mr). (40) 
Jo ^ 

and 

Jo r 

Multiplying equation (39) by ^^MSll which is equal to ^^^^ and integrating gives 

n = n^ + ^j\^ATrr''dr, (42) 

where Q refers to the usual gravitational potential energy of the system. In the end, using equation 
(37) leads to the generalized virial theorem in palatini formalism of f{Tt) gravity 

pi 

2K + n-n^+ / —pA7rr^dr = 0. (43) 
Jo 

Recalling that F = ^ in scalar-tensor representation of the theory, the virial theorem may be written 
in the form 

2iC + J)-0</,+ / j-pAi^r^dr = (44) 
Jo 20 

The fourth term on the left hand side of the virial equation is originated from the relativistic Boltz- 
mann equation in Palatini formalism. This correction term does not exist in metric variational 
approach. To represent the virial theorem in an alternative form, we write equation (43) as 

-R T,/f(„\A^/r rR M^{r)dM F' ^dM 



Jo r Jo r Jq 



2FG r 



= 0, (45) 
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or 



Jo r Jq 



dM 



r 



= 0. (46) 



This equation becomes simpler by the definition M^^"' = — ^qt"^ 

2K-ar^^^-arMr— = o. (47) 

Jo r Jo ^ r 

It is convenient to introduce the radii Ry and 

Jo -^dM{r) 



and 



(M^^"')2 



In addition the virial mass My is defined as [13] 



2if = ^. (50) 



Substituting these definition in equation (47) yields 

My , , {Ml^^fRy 

For most of the observed galactic clusters, the relation My jM > 3 is true. Therefore one can easily 
approximate the last equation 

My (M^n^Ry 

'W ^ Mm^ ■ 

One of the motivations of Palatini fiTZ) gravity models is the possibility of explaining dark 
matter. We noted above that some geometric terms appear in the Einstein field equations which 
could effectively play a role in gravitational energy. These geometric terms may be attributed to a 
geometric mass at the galactic or extra galactic level. They can be interpreted as dark matter which 
originate from modified gravity theory. On the other hand, dark matter provides the main mass 
contribution to clusters. It means that one can ignore the mass contribution of the baryonic mass in 
the clusters and estimate the total mass of the cluster by « M'^^^''" . We also know that the virial 
mass is mainly determined by the geometric mass, so that the geometric mass could be a potential 
candidate for the virial mass discrepancy in clusters. As a result, we conclude that 

M^^^ « My « Mtof (53) 

Therefore, equation (52) can be written as 

R 

My^M-^. (54) 
Ry 

This shows the the virial mass is proportional to the normal baryonic mass of the cluster, whose 
proportionality constant has geometrical origins. 
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5 Astrophysical applications 

5.1 Typical values for the virial mass and radius 

Virial mass density can be written as pv = ^j^g , where My and Ry are the virial mass and 
radius respectively. Astrophysical observations together with cosmological simulations show that the 
virial mass is a measure of a fixed density such as a critical density, pdz) at a special red-shift. 
It means that the virial density can be represented as pv = Spc{z), where 6 ~ 200. As is well 
known, Pc{z) = h'^{z)3HQ /8ttG. The Hubble parameter is therefore normalized to its local value, 
h'^{z) = f^m(l + -z)^ + where ^}„i and Qx ai'c the mass density and dark energy density parameters 
respectively [23]. By knowing the integrated mass of the galaxy cluster as a function of the radius, 
one can estimate the appropriate physical radius for the mass measurement. The commonly used 
radii are r2oo or rsoo- These radii lie within the radii corresponding to the mean gravitational mass 
density of the matter {ptot) = 200pc or 500pc- A useful radius is r2oo to find the virial mass. By 
studying the values of r2oo for various clusters, one can deduce that a typical value for r2oo is 2Mpc 
approximately. The corresponding masses for these radii are defined as M200 and Msoo- In general. 
My = M200 and Ry = ?'2oo are assumed as a virial mass and radius [21]. 



5.2 Geometric meiss and geometric radius from galactic cluster observations 

Intra cluster gas has an important contribution to the baryonic mass of the clusters of galaxies. The 
following equation provides a reasonably good description of the observational date [13, 21] 

2 -— 

Pg = Po (^1 + , (55) 

where Vc is the core radius and pQ and (3 are constants. If we assume that the observed X ray 
emission from the hot, ionized intra cluster gas is in isothermal equilibrium, then the pressure pg of 

the gas satisfies the equation of state pg = {j^^^ Pg, where ks is the Boltzmann constant, Tg is the 

gas temperature, p ~ 0.61 is the mean atomic weight of the particles in the cluster gas and nip is 
the proton mass [13, 21]. Using Jeans equation [1], the total mass distribution can be obtained as 
[12, 13, 21] 

Substitution of the mass density of the cluster gas in equation (55) gives the total mass inside the 
cluster [13, 21] 

'ikBl3Tg 

Now, f{TZ) theory tells us that the total mass of the cluster is Mtotij") = Air j^ipg + P(p)r'^dr, where it 
is assumed that the main part of the baryonic mass of the cluster is in the form of intra cluster gas. 
Therefore, Mtat{r) satisfies the relation 

= ATrr'pgir) + ATrr'p^r). (58) 

Since we have estimated the quantities Mtot{r) and pg, the expression for the geometric mass density 
can be readily obtained 

Urup (ri + r^) ^1 + ^ 
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In the limit r 3> Vc, the geometric density takes the simple form 

47rp<^(r) 

Also, using equation (59), we can easily write the geometric mass as 



finip 



GM,{r) = An r r'p,{r)dr = " - AnGpo f " t a. 

7o ^rnj, i + Lc Jo (1 + ?^)^ 

' c 

In the limit r ^ rc, the geometric mass takes the simple form 



GMJr) 



prrip 3(1 - (3) 



(60) 



(61) 



(62) 



Observations show that the intra cluster gas has a small contribution to the total mass [13, 22, 
21, 23, 24]. This means that the gas density and mass contributions can be neglected compared to 
the geometric density and mass, hence 



Now, it is easy to derive the metric tensor component inside the cluster. Prom equation (36), we 
deduce that r'^v' = -2GMtotir) « -2GM^{r), which leads to 

e'' « ar2«, (65) 
where Ci, is the integration constant and s is defined as 

s = -^-MIl. (66) 

The other coefficient in the metric tensor, e~^, can be estimated approximately as 

_^ ^ 2GMtot 



r 



(67) 



which is the standard expression which, assuming « M^, may be written as e ^ m 1— ^^^"^^ = 
1 + 2s. Of course, there is no guarantee as to the accuracy of this approximation inside the cluster. 
As a result, astrophysical observations can be helpful in deriving the metric components inside the 
cluster in the context of f{TZ) theory . However, how can we estimate the upper bound for the cut 
off of the geometric mass? This happens in a special point where the decaying density profile of the 
geometric density associated with the cluster becomes smaller than the average energy density of the 
universe. We assume that the two densities become equal at the radius R"^. In other words, we 

can write p^(it'^'') = Puniverse- Let us set puniverse = pc = |^ = 4.6975 X 10'^° hl^g/cm'^ where 
H = 50h^okm/Mpc/s [13, 21]. By using equation (63), we obtain 

We may also derive the total geometric mass in equation (64) corresponding to this radius as 

3 

= M,(<)) = 4.83 X lO^V^/^ (^) ' Ko'Me. (69) 

This value is consistent with the mass distribution observations in clusters. Within f{Tl) gravity, 
geometric mass effects are beyond the virial radius, which is about a few Mpc. 
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5.3 Radial velocity dispersion in galactic clusters 

We can write the virial mass in terms of the characteristic velocity dispersion ai as [24] 

Mv = ^<TjRv, (70) 

where Strf = cr^. Let us consider an isotropic velocity dispersion. Then, we have (u^) = (u^) + (uq) + 
{ufp) = 3(n^) = 3(7^ where cr^ is the radial velocity dispersion. Radial velocity dispersion relation for 
clusters in f{TZ) gravity can be obtained from equation (34) as 

-(^i) + -p- + -p^O, (71) 

from which one can deduce a relation for v'. Also equation (36) yields a relation for v' 

1^' = -^[2GM^{r) + 2GM{r) + 2c], (72) 

where c is an integration constant. By eliminating i/' from the last two equations, we obtain 

^(P^r) = -^P+ ^[GM^r] + GM{r) + c]p. (73) 
Integration now gives the solution as follows 

^'^"■^ = ~pl ^''^^ ^~pl ^^^f^'^ + ^^^'^ + + f C^4) 

We can apply equation (74) to a special case where the density p is chosen as [13] 

Pir) = por-^, (75) 

with pq and 7 7^ 1,3 being positive constants. As p is the normal matter density inside the cluster, 
it yields the normal matter mass profile M(r) = 47rpo?'^~'^/(3 — 7). According to equation (64), the 
geometric mass GM<^(r) « qor, where qo = SksPTg/ jimp. The radial velocity dispersion for 7 7^ 1,3 
will be 

(T^{r) = —r^ / — r 'dr — -r ' r'. (76) 

' J 2F 7 (7-l)(3-7) ^ + ir po ^ ^ 

For 7 = 1 we find 

/F' c d 
—r^^dr - go + 2TrGpor Inr - — H r, (77) 
zl< zr po 

and for 7 = 3 we have 

„3 f F'-3r qo , , A Cl C' 3 



<{r) = -r'J dr-^ + 7rGpo[lnr + -j--- + -r-'. (78) 

The observed data can usually be translated to a specific function for velocity dispersion relation. 
Then one can compare the observed velocity dispersion with prediction in modified /(T?.) gravity to 
compare the diflFerent theoretical scenarios. 
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6 The Lagrangian 



We recognized the possibility of finding the metric tensor components in f{TZ) gravity using the virial 
theorem in the previous sections. One of the other consequences is to find the lagrangian f{TZ) of 
the theory, in other words, the form that f{TZ) can take. We start from the field equations in the 
standard representation. Equations (20) and (21) then yield 



2 r 2 F 



Also, from Equations (22) and (23) we have 



1 



3 /F 



/\ 2 



4 V F 



IF" 
+ 27^ + 



u'-X' l\ F 
4 r J f"' 



Equation (2) with /j, = v = 0, using equation (10) leads to 



f = Fe 



-A 



/2 



2u' 



2 + 



3z^'-A' 2\ F F" 
2 '^rjy^^ 



Now, we can write the Ricci scalar from equations (4) and (11) as follows 



(79) 



(80) 



(81) 



1 /F' 



■/\ 2 



- — + — + 



F" 



F' 

+ - I I ^ 



(82) 



Here, R is derived in the context of metric formalism. In fact, we have replaced all the parameters in 
terms of those appearing in the metric formalism which would then lead to the metric Ricci scalar. In 
the derivation of the above equations, we have neglected the contribution from the baryonic matter. 
Equation (65) together with the assumption =const. inside the clusters then become reasonable 
to use. Now, equation (79) can be written as 



FF" - -F'F - -F2 + ^f2 = 0, 



(83) 



for which a solution is 



F = 



4(s2 + 6s + i)r(-'~i-Vs-'+(is+i) 



(84) 



where c is an integration constant and we set c = 1. The component can be derived from Equation 
(80) by using the expression for F 



l + s-s^ + 



(5 -3s- y/s'^ + 6s + 1) . 



Now, equation (81) gives 



/ 



=>-A 



F r2 



-2s - 2s2 + (s + 1 + Vs2 + 6s + l)(-2s + \/s^ + 6s + 1) 



and equation (82) results in 



R 



(s + 1 + + 6s + l)(-5s + 3 + + 6s + 1) , , 2 
^ - 4s - 4s 



(85) 



(86) 



(87) 



Recall that TZ is the Ricci scalar in Palatini formalism whereas R is the Ricci scalar in metric formal- 
ism. Finally, we can easily obtain f as a function of R as follows 



s+3 + V.s-t(>st1 

f = f0R ^ : 



(88) 
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where 



4(s2 + 6s + 1) -2s - 2s2 + (s + 1 + v's2 + 6s + l)(-2s + Vs^ + 6s + 1) 



s+3+V5^+6s+l 



e^^ 2 . (89) 



One may also obtain the action in terms of the Ricci scalar in Palatini formalism. In either case, 
knowing the physical parameters such as gas density and temperature, one can deduce the action for 
modified gravity. 



7 Conclusions 

Virial theorem is a convenient tool with which to derive the mean density of galaxy clusters and it 
can therefore predict the total mass of clusters. We have used the relativistic Boltzmann equation 
within the context of Palatini f{TZ) theory to derive the virial theorem. To write the field equations, 
we have applied the standard cosmological averaging for the energy-momentum tensor. The virial 
mass is mainly determined by the geometric mass which is associated to the geometrical terms in the 
gravitational field equations. The generalized virial mass implicitly includes the effects of dark matter 
and may therefore be used to describe the dynamics of clusters. If equations (54) and (68) accompany 

the assumption « can estimate the virial mass as My ^^■^'^V^i^^^)^ f^5o n^;;^4p^ 

which, by having the physical parameters of a cluster, leads to a specific value. In fact the virial 
mass can be approximated to give the total mass of clusters. In this paper, we have found that the 
virial mass or the total mass of clusters can also be obtained from the velocity dispersion relations. 
However, uncertainties in observational data suppresses the exact value of the virial mass. It seems 
that the gravitational lensing of light in f{Tt) theories gives more exact values for the total mass of 
clusters. In any case, since the virial theorem results can be compared with observational data, one 
can apply it to different theories to test their viability. We also derived, in effect, the Lagrangian of 
the theory in terms of the metric Ricci scalar. Finally, we discussed the limitations of our approach 
which is rooted on the assumption that there exists an averaged metric which is the outcome of 
putting on the right hand side of the field equations an averaged energy momentum tensor. 
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